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ABSTRACT 

We have observed CH absorption lines (7 = 3/2,N =1<—J=l/2,N=1) against the continuum source Sgr B2(M) using the HerschelfHIFl 
instrument. With the high spectral resolution and wide velocity coverage provided by HIFl, 3 1 CH absorption features with different radial veloc- 
ities and line widths are detected and identified. The narrower line width and lower column density clouds show 'spiral arm' cloud characteristics, 
while the absorption component with the broadest line width and highest column density corresponds to the gas from the Sgr B2 envelope. The 
observations show that each 'spiral arm' harbors multiple velocity components, indicating that the clouds are not uniform and that they have in- 
ternal structure. This line-of-sight through almost the entire Galaxy offers unique possibilities to study the basic chemistry of simple molecules in 
diffuse clouds, as a variety of different cloud classes are sampled simultaneously. We find that the linear relationship between CH and H2 column 
densities found at lower Ay by UV observations does not continue into the range of higher visual extinction. There, the curve flattens, which 
probably means that CH is depleted in the denser cores of these clouds. 

Key words. ISM: abundances — ISM: molecules 
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1. Introduction 

While physical and chemical conditions of molecular clouds as- 
sociated with star-forming regions are widely explored by ob- 
serving various molecular species, much less is known about the 
chemistry of the bulk of molecular clouds in the spiral arms 
of our Galaxy. These have mostly been studied only in the 
CO molecule, with some exceptions. These exceptions consist 
of clouds that happen to lie along the line-of-sight to a bright 
Galactic or extragalactic continuum source, and can be observed 
in absorption. Examples for spiral arm clouds probed by gas ab- 
sorption are the Sgr B2, W49, W51, and Cas A millimeter con- 
tinuum sources, where the chemistry was studied mostly through 
low-energy lines of molecular species (e.g.. Greaves & Williams 
1994; Greaves & Nyman 1996; Gerin et al. 2010; Menten et al. 
2010; Lis et al. 2010; Schilke et al. 2010; Tieftrunk et al. 1994; 
Wyrowski et al. 2010). Liszt and Lucas (2002) also made impor- 
tant contributions to this field through observations of absorption 
line clouds against quasars. The observations show that most of 
the absorption signal arises from diffuse/translucent clouds with 
low gas density and excitation temperature, suggesting that the 
ambient ultraviolet field plays a prominent role in determining 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



the chemistry (Lis et al. 2010; Gerin et al. 2010; Neufeld et al. 
2010). 

Simple linear molecules are fundamental constituents when 
studying interstellar chemistry, since they are building blocks 
for more complex species (e.g., Liszt 2009). The very light hy- 
dride radical CH (methylidyne) has been observed in diffuse 
clouds in front of bright OB stars at optical wavelengths and 
in dark clouds at radio wavelengths, showing that CH column 
density is linearly correlated with the optical extinction (A,,) for 
A,, less than 4'" (e.g., Lang & Wilson 1978; Federman 1982; 
Sheffer et al. 2008). Observations and chemical models also in- 
dicate that CH abundance is lower by a factor of 100-1000 in 
dense clouds and is mainly located in the surface layers of dense 
clouds (Viala, Bel & Clavelet 1979; de Jong, Genzel et al. 1979; 
Boland & Dalgarno 1980; Goicoechea, Rodriguez-Fernandez & 
Cernicharo 2004; Polehampton et al. 2007) or low column den- 
sity clouds (Mattila 1986; Magnani, Lugo & Dame 2005). 

The prominent Sagittarius B2 star-forming region is very 
close (~ 130 pc) to the Galactic center and at a distance of ~8 
kpc from the Sun (Reid et al. 2009). The dense cores contained 
in it, Sgr B2(N) and Sgr B2(M), are well-studied massive star- 
forming regions. Sgr B2(M) is stronger at submm wavelengths 
and also less contaminated by emission lines than its neigh- 
bor Sgr B2(N) (Goldsmith et al. 1990; NummeHn et al. 1998), 
making it the target of choice for absorption studies. HIFI, the 
Heterodyne Instrument for the Far-Infrared (de Graauw et al. 
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2010) onboard the Herschel Space Observatory (Pilbratt et al. 
2010) is an ideal instrument for making such observations. In 
this Letter, we present high spectral-resolution HlFl observa- 
tions of CH towards Sgr B2(M). 



2. Observations 

A full spectral scan of Sgr B2(M) in the HIFl band la, cov- 
ering a frequency range from 479.6 to 560.3 GHz was car- 
ried out on March 1, 2010, and was pointed at q'(J2000) - 
17H7"'20.2' and 5(J2000) =-28°23'05.0", using the dual beam 
switch (DBS) mode as part of the guaranteed time key program 
HEXOS: Herschel/HIFI Observations of Extraordinary sources: 
The Orion and Sagittarius B2 Star forming Regions (Bergin et 
al. 2010). The DBS reference beams lie approximately 3' east 
and west (i.e. perpendicular to the roughly north-south elonga- 
tion of Sgr B2). The 1.1 MHz spectral resolution of the Wide 
Band Spectrometer (WBS) corresponds to a velocity resolution 
of ~0.6 km s ' over a 4 GHz IF bandwidth. The data have been 
processed through the standard pipeline released with version 
2.9 of HIPE (Ott et al. 2010), including removal of the standing 
waves, and subsequently exported to FITS format. The double- 
sideband (DSB) spectra were deconvolved into single-sideband 
spectra, including the continuum (Comito & Schilke 2002) us- 
ing the IRAM GILDAS package. The systematic velocity of the 
source was corrected to the local stand of rest (LSR) frame. Both 
H and V polarization data were obtained. The data presented 
here are averages of the two polarizations with equal weighting. 
The HlFl beam size at 532 GHz is ~37", with a main beam effi- 
ciency of 0.69. 



3. Results 

The low-frequency radio transitions of CH are between A-type 
doublet states, while spin-rotational as well as rotational transi- 
tions fall in the far-infrared and submillimeter bands (Brown & 
Evenson 1983). Early observations mainly focused on the low- 
frequency radio and far-infrared lines, the latter with relative 
low spectral resolution. The lowest spin-rotational transitions of 
CH at 532.8 and 536.8 GHz were measured in the laboratory 
by Amano (2000), thus lying in a range of strong atmospheric 
water absorption, which makes their observation with ground- 
based telescope impossible. The level populations for the 532.8 
and 536.8 GHz transitions are expected to have a low excitation 
temperature in most clouds, with the possible exception of the 
Sgr B2 envelope, and thus provide a good measure of the total 
column density from absorption measurement. 

Six hyperfine components of CH in the ground electronic 
state {X ^H), J = 3/2,A^ = 1 ^ / = 1/2,A? = 1 near 
532.8 (1-1- ^ 1-,2h- ^ 1-, and \+ <r- 0-), and 536.8 GHz 
(2- <— l-H, 1- <— l-H, and 1- «— On-) are included in the 
HIFl band la observations of Sgr B2(M). The CH spectra near 
532.8 and 536.8 GHz obtained with HlFl are presented in Fig. 
1, with multiple velocity components seen in absorption. The 
observations of CH at 7 = 3/2, A^ = 2<-7=l/2,A^ = 2 
in far-infrared (~ 2010 GHz) with a velocity resolution of >30 
km s"' towards Sgr B2 using the Infrared Space Observatory 
(ISO) and the Kuiper Airborne Observatory (KAO) 91.4 cm tele- 
scope showed absorption centered on 65 km s ' with an unusu- 
ally large line width of -200 km s"' (Stacey, Lugten & Genzel 
1987; Goicoechea, Rodriguez-Fernandez & Cernicharo 2004; 
Polehampton et al. 2007), which samples the gas from both the 
Galactic center and spiral arm clouds along the light of sight. 



The HIFl CH spectra are spectrally resolved into multiple ve- 
locity components with different line intensities and widths. The 
strong CH absorption at ~64 km s ' with the largest line width 
clearly corresponds to the envelope of the Sgr B2(M) cloud. In 
addition to the Sgr B2(M) absorption cloud, the radial velocities 
of the other absorption components are consistent with those of 
the spiral arm clouds and kinematic features associated with the 
Galactic center (e.g.. Greaves & Williams 1994; Menten et al. 
2010). HIFl observations of H20^ and H2O in the spiral arm 
clouds against the Sgr B2(M) continuum show somewhat dif- 
ferent spectral profiles than CH (Schilke et al. 2010; Lis et al. 
2010). Probably each molecular species samples slightly differ- 
ent regions within various clouds along the line of sight. 
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Fig. 1. CH spectra with velocity scale referring to the hyper- 
fine fines at 536.761 145 and 532.723926 GHz. The upper panel 
shows the spectrum at 536.8 GHz, the lower panel shows the 
spectrum at 532.8 GHz. The green curve in the bottom panel 
represents the total opacity as a function of LSR velocity, ob- 
tained by linearly adding up the opacities of the hyperfine com- 
ponents. In each panel, the black curve is the observed spectrum, 
the LTE fit is shown as the red curve, and the blue curve shows 
the hyperfine pattern for the transition. 



Absorption-line column densities can be calculated from 
the observed line-to-continuum ratio by Gaussian fitting to the 
lines, if the absorption spectrum comes from a single molecular- 
species without hyperfine components, as long as the spectrum 
is not contaminated by emission from other species. We have 
not identified any strong interfering lines in the velocity range of 
the 532.8 and 536.8 GHz CH absorption. However, the hyperfine 
structure of CH makes it difficult to uniquely assign the absorp- 
tion to distinct velocity components and calculate their column 
densities. Taking the hyperfine components and contamination 
from other species into account, and using the LTE approxima- 
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Table 1. Model-fitting results and abundances 
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a, The H2 column densities are decoded from HCO* column densities 
(Menten et al. 2010). ft, SgrB2(M); c, Scutum arm; d, Sagittarius arm; 
e, Galactic center; / Norma arm (3 and 4 kpc arms). 



tion, we modeled the CH spectra by employing XCLASSQ and 
the automated fitting routine provided by MAGIXq We adopted 
2.73 K as rotation temperature for all absorption components, 
which is questionable only for the Sgr B2 component itself, since 
the envelope is hot and dense enough to excite CH. A reliable es- 
timate of the CH column density and abundance in this compo- 
nent has to await detailed radiative transfer modeling. Keeping 
the velocity offsets fixed, we used 31 absorption components 
and 1 emission component, as for C-C3H2 and H'^CO^ cases 
(Menten et al. 2010) for the CH model fitting. We assume that 
the opacity at any velocity is a sum of contributions from all CH 
hyperfine components and that the absorbing layers cover the 
continuum completely. The best fits to the spectra along with the 
line opacity as a function of velocity are shown in Fig. 1. The 
column densities and line widths of the 3 1 absorption compo- 
nents of the best-fit model are given in Table 1 . Based on the 
results of Menten et al. (2010) and Greaves & Williams (1994), 
the corresponding name of the spiral arm clouds for each veloc- 
ity component is also given in Table 1 . 



' We made use of the myXCLASS program (https://www.astro.uni- 
koeln.de/projects/schilke/XCLASS), which accesses the CDMS 
(Muller at al. 2001, MuUer et al. 2005; http://www.cdms.del and JPL 
(Pickett et al. 1998; http://spec.jpl.nasa.gov I molecular data bases. 

2 https://www.astro.uni-koeln.de/projects/schilke/MAGIX 



The column densities range from 2.3x10 to 
5.2x10"* cm-^ The Sgr B2 cloud (-64 km s-') has the 
largest line width and highest column density. Most of the com- 
ponents have column densities of a few times 10'^ cm-^. The 
mean column density is 5.4x10'-' cm-^, excluding Sgr B2(M). 
This is similar to the Herschel/HIFI observations of CH towards 
the massive star-forming region NGC6334 (Van der Wiel et al. 
2010), in which the NGC6334 envelope has the highest column 
density of l.TxlO'"* cm~^, while foreground clouds have a 
column density of several 10'^ cm-^. The CH column densities 
in our observations are similar to those in diffuse clouds and 
dark nebulae (e.g., Federman 1982; Liszt & Lucas 2002; 
Mattila 1986). Each spiral arm cloud harbors several velocity 
components with different line widths and column densities. 
The velocity structure in each spiral arm cloud suggests that the 
clouds are not uniform and have internal structure. 



4. Discussion 

For absorption lines, both the observed antenna temperature and 
column density of a specific molecule are related to optical depth 
and excitation temperature. The optical depths in this work are 
derived from the observed line-to-continuum ratio, therefore the 
uncertainties due to the absolute flux calibration will not affect 
the values of the optical depth. The brightness temperature of the 
background continuum is determined from the baseline offset, 
which is fairly reliable for HIFI (Bergin et al. 2010). The only 
assumption is that the absorption completely covers the contin- 
uum. The spiral-arm cloud sizes are approximately 3 pc corre- 
sponding to 150" at a distance of 4 kpc (Greaves & Williams 
1994), which is more than the size of the Sgr B2(M) continuum. 
The same set of parameters in the LTE model fits both the spec- 
tral features near 532.8 and 536.8 GHz well. The typical optical 
depths are 0.5, while the largest optical depth is ~1 at 64 km s-' 
(see Figure 1). Thus, our LTE fitting gives a reliable assessment 
of the column densities. 

Determining of the fractional abundance of a molecule rel- 
ative to H2 requires the column densities of both the molecules 
in question and H2. A direct measurement of H2 column den- 
sity is difficult. HCO^ has been shown to have an approxi- 
mately constant abundance relative to H2 in a wide range of en- 
vironments, including diffuse clouds, low- and high-mass star- 
forming regions, and compact extragalactic continuum sources 
(e.g.. Greaves & Nyman 1996; Lucas & Liszt 1996; van 
Dishoeck et al. 1993); therefore, HCO^ can be taken as a proxy 
of H2. However, HCO^ itself is optically thick in most veloc- 
ity components. Following Menten et al. (2010), we use the op- 
tically thin H'^'CO^ column densities from Table 4 of Menten 
et al. (2010) to determine the H2 column density and fractional 
abundance of CH relative to H2 for each velocity component. 
This procedures takes measured '^C/'-'C ratios in the various 
spiral-arm clouds into account and assumes a constant HCO^/H2 
ratio of 5x 10-''. The fractional abundances of CH relative to H2 
are listed in the fifth column of Table 2. The derived H2 column 
densities range from 1x10^' to 3x10^^ cm -^ with a total col- 
umn density of 1.6x10^'' cm -^ except for Sgr B2(M), which are 
roughly consistent with the values derived from '^CO in spiral- 
arm clouds (Greaves & Nyman 1996; Neufeld et al. 2000). Most 
of the gas components have a CH fractional abundance of a few 
times 10-*^ with an average abundance of 2.1x10"^, excluding 
the Sgr B2(M) cloud with a value of 9.6x10-'" because the ex- 
citation there is not straightforward. 

Figure 2 presents the relationship between the CH and H2 
column densities. A least-squares fit (excluding the Sgr B2(M) 
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Fig. 2. Correlation of column densities of CH and H2 in a log- 
log plot. H2 column densities are derived from H^^CO^ column 
densities (Menten et al. 2010). The least-squares fit to the data is 
indicated by the red line. The green line shows the slope of 0.97 
fitted to the diffuse clouds in a range of 19.2 < logA^(H2) < 21.15 
(Sheffer et al. 2008). 



envelope point) gives a slope of a = 0.38 + 0.07 and an in- 
tercept of b = 5.81 + 1.64 (logA^(CH) = alogNiUi) + b) 
with a correlation coefficient of 0.8. The fits to diffuse clouds 
(19< logA^(H2) <21.15) give a slope of ~1 with a correlation 
coefficient of 0.99 (Federman 1982; Sheffer et al. 2008), sug- 
gesting that A^(CH) correlates linearly with A^(H2). In Figure 
2, most of the velocity components in our observations have 
higher H2 column densities compared to those of diffuse clouds 
(Sheffer et al. 2008), and the linear relationship breaks down 
with N(H2) > 10^' cm"', i.e. at visual extinctions beyond 
1"', where has a flatter slope of 0.38. The observations of CH 
and H2CO towards a few dense molecular clouds associated 
with massive star formation regions (Genzel et al. 1979) have 
shown that A^(CH)/A^(H2CO) decreases with increasing H2CO 
column density, suggesting that CH abundance is lower in dense 
clouds. The relatively flat slope at higher H2 column densities 
and corresponding lower fractional abundances in our observa- 
tions corraborate the evidence that CH abundance is lower in the 
inner part of the clouds and predominantly traces lower density 
environments. As a word of caution, we note that this hinges on 
the assumption of constant abundance of HCO^ and has to be 
tested with other tracers for H2 column densities. 

Previous observations have also shown that the CH column 
density in photodissociation regions (PDR) is not consistent with 
the linear relationship (Sheffer et al. 2008). Our results appear 
to agree with the predictions of both turbulent dissipation re- 
gions (TDR) and PDR models (see Fig. 8 of Godard, Falgarone 
& Pineau des Forets 2009) with an H2 volume density less than 
200 cm"^. The chemical network suggests that CH, CH^, CH2, 
C, CH3, CHj, C+, C2+, etc. are tightly related by chemistry 
(Godard, Falgarone & Pineau des Forets 2009). Comparing the 
column densities of all the hydrides and molecular ions in the 
absorption line clouds, as provided by HIFI, will enable us to 
investigate the relative contributions of photon dominated and 
turbulence dominated chemistry, and distinguish PDR and TDR 
models. 



References 

Amano, T., 2000, ApJ, 531, L61 

Bergin, E. A., Phillips, T. G., Comito, C, et al., 2010, this volume 

Brown, J. M., & Evenson, K. M., 1983, ApJ, 268, L51 

Comito, C, & Schilke, P, 2002, A&A, 395, 357 

de Graauw, Th. et al. 2010, A&A, in press (Herschel special issue) 

de Jong, T, Boland, W., & Dalgarno, A., 1980, A&A, 91, 68 

Federman, S. R., 1982, ApJ, 257, 125 

Genzel, R., Downes, D., Pauls, T., Willson, R. F, Bieging, J., 1979. A&A, 73, 

253 

Gerin, M., De Luca, M., Black, J., et al. 2010, A&A, this volume 

Goldsmith. PF, Lis, D. C, Hills, R., Lasenby, J., 1990, ApJ, 350, 186 

Godard, B., Falgarone, E., & Pineau Des Forets, G., 2009, A&A, 495, 847 

Goicoechea, J. R., Rodriguez-Fernandez, N. J., & Cemicharo, Jose, 2004, ApJ, 

600,214 

Greaves, J. S., & Williams, R G., 1994, A&A, 290, 259 

Greaves, J. S., & Nyman, L.-A., 1996, A&A, 305, 950 

Lang, K. R., & Wilson, R. F, 1978, ApJ, 224, 125 

Leung, C. M., Herbst, E., & Huebner, W. F, 1984, ApJS, 56, 231 

Lis, D. C, Phillips, T. G., Goldsmith, P F, et al., et al. 2010, this volume 

Liszt, H. S., & Lucas, R., 2002, A&A, 391, 693 

Liszt, H. S., 2009, A&A, 508, 783 

Lucas, R., & Liszt, H., 1996, A&A, 307, 237 

Magnani, L., Lugo, S., & Dame, T M., 2005, AJ, 130, 2725 

Mattila, K., 1986. A&A, 160, 157 

Menten, K. M., Wyrowski, F, Belloche, A., et al., 2010, submitted to A&A 

Miiller, H. S. P, Thorwirth, S., Roth, D. A., & Winnewisser, G., 2001, A&A, 

370, L49 

Muller, H. S. P, Schloder, F, Stutzki, J., & Winnewisser, G., 2005, J. Mol. 

Struct, 742, 215 

Neufeld, D. A., Ashby, M. L. N., Bergin, E. A. et al., 2000, ApJ, 539, LI 11 

Neufeld, D.A., Sonnentrucker, P., Phillips, T. G., et al., 2010, this volume 

Nummelin, A., Bergman, P, Hjalmai'son, A., et al., 1998, ApJS, 117, 427 

Ott, S., et al., 2010, in ASP Conference Series, Astronomical Data Analysis 

Software and Systems XIX, Y. Mizumoto, K.-I. Morita, and M. Ohishi, eds., in 

press 

Pickett, H. M., Poynter, R. L., Cohen, E. A., et al., 1998, J. Quant. Spectrosc. 

Radiat. Transfer, 60, 883 

Pilbratt, G., Riedinger, J. R., Passvogel, T., et al., 2010, A&A, in press 

(Herschel special issue) 

Polehampton, E. T, Baluteau. J., Swinyard. B.M., et al., 2007, MNRAS, 377, 

1122 

Reid, M. J., Menten, K. M., Zheng, X. W., Brunthaler, A., Xu, Y, 2009, ApJ, 

705, 1548 

Schilke, P, Comito, C, Muller, H. S. P, et al., 2010, this volume 

Sheffer, Y, Rogers, M., Federman, S. R., et al., 2008, ApJ, 687, 1075 

Stacey, G. J., Lugten, J. B., & Genzel, R., 1987, ApJ, 313, 859 

Tieftrank, A., Pineau des Forets, G., Schilke, P, & Wahnsley, C. M., 1994, 

A&A, 289, 579 

Van der Wiel, M. H. D., van der Tak, F F. S, Lis, D. C, et al., 2010, this volume 

van Dishoeck, E. F, Blake, G. A., Draine, B. T., Lunine, J. I., 1993, in 

Protostai's and planets 111, ed. E. H. Levy & J. 1. Lunine (Tucson, AZ: Univ. 

Arizona Press), 163-242 

Viala, Y P, Bel, N., & Clavel, J., 1979, A&A, 73, 174 

Wyrowski, F, Menten, K. M., Guesten, R., Belloche, A, 2010, A&A, in press 



S.-L. Qin et al.: Spiral arm clouds by CH absorption lines 



I. Physikalisches Institut, Universitat zu Koln, Ziilpicher Stn 77, 

50937 Koln, Germany 

e-mail: qin@phl.uni-koeln.de 

Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, 

53121 Bonn, Germany 

California Institute of Technology, Cahill Center for Astronomy and 

Astrophysics 301-17, Pasadena, CA 91125 USA 

Department of Astronomy, University of Michigan, 500 Church 

Street, Ann Arbor, MI 48109, USA 

LERMA & UMR8112 du CNRS, Observatoire de Paris, 61, Av. de 

rObservatoire, 75014 Paris, France 

Centre d'etude Spatiale des Rayonnements, Universite de Toulouse 

[UPS], 31062 Toulouse Cedex 9, France 

CNRS/INSU, UMR 5187, 9 avenue du Colonel Roche, 31028 

Toulouse Cedex 4, France 

Laboratoire d'Astrophysique de I'Observatoire de Grenoble, BP 53, 

38041 Grenoble, Cedex 9, France. 

Centro de Astrobiologia (CSIC/INTA), Laboratiorio de Astrofisica 

Molecular, Ctra. de Torrejon a Ajalvir, km 4 28850, Torrejon de 

Ardoz, Madrid, Spain 

LERMA, CNRS UMR8112, Observatoire de Paris and Ecole 

Normale Superieure, 24 Rue Lhomond, 75231 Paris Cedex 05, 

France 

LPMAA, UMR7092, Universite Pierre et Marie Curie, Paris, France 

LUTH, UMR8102, Observatoire de Paris, Meudon, France 

Jet Propulsion Laboratory, Caltech, Pasadena, CA 91 109, USA 

Departments of Physics, Astronomy and Chemistry, Ohio State 

University, Columbus, OH 43210, USA 

National Research Council Canada, Herzberg Institute of 

Astrophysics, 5071 West Saanich Road, Victoria, BC V9E 2E7, 

Canada 

Infrared Processing and Analysis Center, California Institute of 

Technology, MS 100-22, Pasadena, CA 91 125 

Canadian Institute for Theoretical Astrophysics, University of 

Toronto, 60 St George St, Toronto, ON M5S 3H8, Canada 

Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, 

Cambridge MA 02138, USA 

National University of Ireland Maynooth. Ireland 

Department of Physics and Astronomy, Johns Hopkins University, 

3400 North Charles Street, Baltimor e, MD 21218, USA 

SRON Netherlands Institute for Space Research, PO Box 800, 9700 

AV, Groningen, The Netherlands 

Department of Physics and Astronomy, University of Calgary, 2500 

University Drive NW, Calgary, AB T2N 1N4, Canada 

MPI fiir Sonnensystemforschung, D 37191 Katlenburg-Lindau, 

Germany 

Observatorio Astronomico Nacional (IGN), Centro Astronomico de 

Yebes, Apartado 148. 19080 Guadalajara, Spain 

Microwave Laboratory, ETH Zurich, 8092 Zurich, Switzerland 

Chalmers University of Technology, SE-412 96 Gteborg, Sweden, 

Sweden/ SRON Netherlands Institute for Space Research, 

Landleven 12, 9747 AD Groningen, the Netherlands 

Department of Astronomy, Stockholm University, SE-106 91 

Stockholm, Sweden 



